clc
clear

theta = deg2rad(30);  
A = [cos(theta) -sin(theta) 0; sin(theta) cos(theta) 0; 0 0 1];
fprintf("Rotation Matrix:\n");
disp(A);

P = [0 2 0]';
fprintf("P refer to A Axes\n");
disp(A*P);
